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Comment on "First order amorphous- amor- 
phous transformation in silica" 

In a recent letter , Lacks presents evidence of a first 
order amorphous-amorphous transition in silica at T = 
0. He bases his conclusions on the calculation of the 
free energy along a path of compression and successive 
decompression of a sample of 108 Si02 units. The free 
energy of the two branches cross each other, and this is 
interpreted as evidence of a first order transition, which is 
suggested to be hidden in the actual simulations (and in 
experiments) due to dynamic arrest. We show that this 
conclusion does not follow from the shown data, since 
qualitatively the same phenomenology is obtained in a 
model where a first order transition does not exist. 
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FIG. 1. (a) Free energy g(P) and volume v(P) forms for 
a ID interaction potential with two equilibrium distances 
(sketched in the inset), (b) The result of averaging over a 
disordered sample. 

The model we use has been applied to study the 
anomalous properties of tetrahedrally coordinated liquids 
of which silica is an example. It consist of particles 
interacting through a spherical, purely repulsive poten- 
tial U{r) as qualitatively depicted in the inset of Fig. 1. 
Its main property is the existence of two different equilib- 
rium distances depending on pressure. We will work out 
here the one-dimensional (ID) case, that can be solved 
exactly, but we have performed numerical simulations in 
more realistic two- and three-dimensional systems, with 
the same qualitative results [|| (see also 1^, Fig. 7). 
For the ID system at zero temperature, the qualitative 


evolution of Gibbs free energy and volume vs. external 
pressure for a pair of neighbor particles is given in Fig. 
1(a). There is a bistability upon compression and decom- 
pression, as indicated in the figure. When an assembly of 
particles is considered, we have to sum over all of them to 
get the result for the whole system. In ID, polydispersity 
should be introduced to avoid 'crystalHzation'. The free 
energy has to be calculated as G{P) = J2i9^{P)j where 
i labels the particles, and g^{P) is the free energy of each 
particle, which is not the same throughout the system 
due to polydispersity (in real glasses, the abundance of 
different local environments plays the role of the poly- 
dispersity we introduce in the ID case). Depending on 
what we are interested in, the compression, decompres- 
sion, or thermodynamic forms oi g^{P) have to be used. 
An analog expression holds for V{P). 

The qualitative final forms of G{P) and V{P) are 
shown in Fig. 1(b). The free energy of compression 
and decompression branches cross at some point (which 
is reminiscent of the microscopic g{P) form) as in Fig. 
2(a) of [Q. However, the thermodynamic G{P) is always 
lower, and smooth for all P, without signs of a first order 
transition. Then we conclude that the evidence presented 
in [|l| for silica is not neccesarely related to a first order 
transition. The definite answer must be sought in sim- 
ulations in which the system is annealed at different P, 
to obtain the thermodynamic form of G{P) and V{P). 
For T > we expect hysteresis in V{P) and G{P) to dis- 
appear in simulations (or experiments) performed at T 
higher than the glass temperature Tg, at which dynamics 
freezes. 

Finishing, we notice that the presence of an attractive 
part in the interaction potential may produce (if strong 
enough) a first order amorphous-amorphous transition in 
the model In this case, numerical simulations at zero 
temperature in a two dimensional system |^ show that 
our model reproduces qualitatively the V{P) behavior 
observed in water [Q. This points out that a first order 
amorphous-amorphous transition-if it cxists-can be di- 
rectly observed in numerical simulations even well below 
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